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ABSTRACT 


Fatigue data generated by three combined bending- torsion fatigue 
reliability research machines at The University of Arizona are proba- 
bilistic-graphically and phenomenologically analyzed. Distributions 
that are applicable to fatigue life and static strength data are dis- 
cussed. Phenomenological justifications for the use of these distri- 
butions are presented. 

It is found that the normal distribution represents the cycles- 
to- failure data at the highest stress levels best. The lognormal 
distribution appears to fit the cycles-to-failure data at the lower 
stress levels best and quite well at all stress levels including the 
highest. 

A regression analysis and least-squares goodness-of-fit test 
was performed for normal and lognormal plots. In most cases, the 
correlation coefficient gave a better fit to the data using the nor- 
mal distribution, but the difference between the two was so slight 
that positive discrimination could not be made. 

From the probabilistic-graphical analysis and the phenomeno- 
logical reasoning, it was concluded that the lognormal distribution 
gave a very satisfactory fit to the cycles-to-failure data at all 
stress levels, that the normal distribution could be used to repre- 
sent the cycles-to-failure at the highest stress levels without any 
loss of accuracy. 

The normal distribution is found to describe the static strength 
distributions best. 

The Weibull distribution was also studied and the probabilistic- 
graphical plots were found not to lend themselves to as good a 
straight line fit 'to the cycles-to-failure data as desired. The plots 
had many. kinks that could not be straightened by adjusting the loca- 
tion parameter. Phenomenologically it was found that the Weibull 
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would be the best distribution if the cycles-to-failure data were 
those of only the failed items in a large sample or of field failures 
because these would be the failures of the weakest in a sample thus 
conforming to the extreme-value Weibull distribution theory. 
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I SUMMARY . 

Three combined-stress fatigue reliability research machines 
have been built and calibrated to test relatively Xsrge specimens 
under combined reversed-bending and steady-torque loading condi- 
tions. These machines are generating. data to be used in deter- 
mining statistical strength surfaces (three-dimensional Goodman 
diagrams) so that specified reliabilities may be designed into 
components subjected to such combined loading using the design 
by reliability methodology. 

This report is based on the results from 170 test specimens 
which yielded cycles-to- failure, stress-to-failure, and endurance 
strength data at various stress levels and at the alternating 
stress to mean stress ratios of «>, 0.70, 0.90, and 0. The speci- 
mens were made of SAE 4340 Steel of Rockwell 35 to 40 hardness 
on Scale C, and processed according to MIS-S-5000B, MIL-H-6875, 
and MIL- 1-6868. *: 

Specimens were tested at five levels of alternating stress 
for the stress ratio of <», i.e., 154,000 psi, 121,000 psi, 

104,500. psi, 87,000 psi, and 78,000 psi. Eighteen specimens were 
tested at each of the four lower stress levels and twelve specimens 
at the highest level. Specimens were tested at four levels of al- 
ternating stress for the stress ratio of 0.70, i.e., 110,500 psi, 
97,500 psi, 76,000 psi, and 70,000 psi. Eighteen specimens were 
tested at each of the three lower levels and twelve specimens at 
the highest level. 

Ten notched and ten unnotched specimens were statically tested 
to failure for the stress ratio of 0, and the ultimate strength 
was determined. 

The Weibull, normal and lognormal distributions were fitted 
to the cycles-to- failure and ultimate strength data. Graphical- 



2 


probabilistic and phenomenological analysis was made of the data 
to decide which statistical distribution best • represents them. 

The normal and lognormal distributions gage good fits to the 
cycles- to-failure data for the stress ratios of «> and 0.70. The 
Weibull distribution, although very versatile, is an extreme-value 
distribution and does not exactly reflect the results of this re- 
search. Phenomenological reasons favored the lognormal over the 
normal distribution to best represent the cycles-to-failure data 
at all stress levels. At the highest stress levels , however, the 
normal distribution can be used to approximate the lognormal dis- 
tribution. This can be justified both probabilistic-graphically 
and phenomenologically. 

The static ultimate strength data of the notched specimens 
.(stress ratio of 0) were found to be best represented by the normal 
distribution both graphically and phenomenologically. 
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II INTRODUCTION 

When a specimen is subjected to an alternating stress, even a 
stress below the static fracture strength, cracks will form and 
propagate to cause rupture or failure. This phenomenon is called 
fatigue, and the rupture is referred to as a fatigue failure. If 
the same level of alternating stress is applied to several speci- 
mens , the scatter in the number of cycles necessary to produce 
fatigue failure is quite large. 

This scatter exceeds experimental error, and many testing pro- 
grams have been executed which show the scatter cannot be explained 
by differing surface finishes , heat treatments , or inhomogeneity of 
the material. Physicists, metallurgists, and applied mathematicians 
have proposed various theories which place a significance on the 
scatter from a statistical viewpoint. 

The fact that this scatter in cycles-to-failure at a constant 
stress level exists has been known for a long time. It has only 
been quite recently, however, that the statistical nature of fatigue 
has been recognized. It is now considered a fundamental and essen- 
tial characteristic of fatigue analysis . 

For three years, fatigue testing, under the combined effects of 
alternating bending and mean torque, has been carried on at The - 
University of Arizona. This testing was done under the direction of 
Dr. Dimitri Kececioglu and the sponsorship of the National Aero- 
nautics and Space Administration. 

The testing has been carried out on combined bending-torsion, 
fatigue reliability research machines designed and built at The 
University of Arizona. These machines are capable of applying and 
maintaining an alternating bending stress and a mean shear stress 
at different levels in a rotating round specimen. The machines and 
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specimens -"are described in detail in previous reports (l, pp. 193- 
257 ) y _(2) . „The test specimens were notched with a theoretical 
stress concentration of 1.45 and were of SAE 4340 Steel. 

__ The . tests were conducted at various alternating bending stress 
levels while holding constant alternating-stress-to-mean stress 
ratios... The. purpose was to determine statistically the effects of 
superposing. steady torque onto bending on the S-N diagram for such 
constant . stress ratios. Testing has been completed for three stress 
ratios: 0.70 (0.90 for endurance) and 0. The data and results 

are given in Tables 1, 2 and 3 and Figures 1 and 2. 

The data can be used to generate statistical distributions of 
cycles -to-failure and a statistical S-N diagram at each stress ratio 
as shown in Figures 1 and 2. This information can also be used to 
determine the strength of the specimens at specific cycles of life 
as shown in Figure 3. 

^ After data for several stress ratios have been gathered, a 
statistical Goodman surface for the fatigue strength of a specimen 
at specific cycles of life for various stress ratios can be gener- 
ated as shown in Figure 4. 

The stress ratio, r , is defined as 

s 


s 



For The University of Arizona research program the alternating stress 

is a bending stress , and the mean stress is a shear stress . Using 

the von Mises-Hencky theory of failure s = a and s - S3 r (2. 

a a m m 9 

p. 87), the stress ratio then becomes 


a 
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A. Stress Ratio of 


00 


To test a ratio of 00 is to test the specimens with pure bending 
stress only and zero shear stress. This was done at five levels of 
alternating stress and the cycles-to-failure were determined. Test- 
ing was also done, using the staircase method (3, pp. 113-114), to 
find the endurance strength of the notched specimens . The results 
are presented and discussed in a previous report (2). 

The alternating stress levels at r = 00 were 

s 

& 

s = 154,000 psi, 

3 1 - • 

s = 121,500 psi, 
a 2 

s = 104,500 psi, ' . 
a 3 

s = 87 ,000 psi, 
a 4 

S a^ = 78,000 psi. 

The cycles-to-failure at these various alternating stress levels are 
given in Table 1. The results are given in Figure 1. 


These stresses are rounded out to the nearest 500 psi. 



B. Stress Ratio of 0.70 . 

To test at a stress ratio of 0.70, various levels of alternating 
bending stress were used with the mean shear stress adjusted to main- 
tain the constant ratio. Four levels of alternating stress were used 
They were 

s = 110,500 psi, 
a l 


s = 97,500 psi, 

s = 76,500 psi, 

a 3 

s = 70,000 psi/' 
a 4 

Using the von Mises-Hencky theory, the shear stress at each level was 
= 88,500 psi, = 80,500 psi, x^ = 65,000 psi, x^ = 57,000 psi, 
to maintain a stress ratio of 0.-70. The cycles-to-failure data are 
given in Table 1. The staircase method was used to determine the 
distribution of the endurance strength at a stress ratio of 0.90. 

The results are given in Figure 2. 



C. Stress Ratio of 0 j 

To determine the distribution for stress ratio of 0, twenty spe- 
cimens were tested to failure under tensile loading. Ten of the 
specimens were regular notched specimens and ten were unnotched with 
a diameter equal to the diameter at the base of the notch. The 
results are given in Tables 2 and 3. 

The strength distribution for a stress ratio of 0 was taken to 
be the ultimate strength distribution for the notched specimens. 

The ultimate is the end of the load carrying ability of the specimen, 
and this is of interest to the design engineer. The ultimate 
strength, as the intersection of the modified Goodman line and the 
abscissa, agrees with the literature (2, pp. 20-23),$+, p. 180), (5, 
p. lTsJjCOjP. 270). Figure 4 reflects the use of this conclusion. 
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_ III. THE DENSITY FUNCTIONS, MOMENTS AND PARAMETERS OF 
STATISTICAL DISTRIBUTIONS 

A. Density Functions 

The cumulative distribution, F(x^), of a random variable x is 
the probability that x assumes a value no greater than some speci- 
fied x. , or 
1 


:F-(x^) = -Probability (x - x^) . 


( 3 ) 


The probability density function f(x) is defined as 


Probability (x. 1x5 x.+Ax.) 

f (x) = lim -±- i — 

Ax^-K) 4x i 


or 


f(x) = ^ CF(x)] ' 


( 4 ) 


It follows from the definition of the probability density function 
that 


f(x)dx = F(x) 


= F(°°) - F( -“) = 1-0 = 1 , 


--00/ 
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or the area under the probability density function is equal to one. 
Also it should be noted that the cumulative distribution can be 
obtained from the probability density function as follows: 


F(t) = 


The probability density function is referred to as the pdf, and the 
cumulative distribution function as the cdf. » - 
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B. Expected Values and Moments of Distributions 

The expected value of a distribution, or the mean of a random 
variable x, is of interest in analyzing distributions. This is 
given by 


'E(x) = 


xf (x)dx = y 


( 6 ) 


The moment of the distribution about the mean, y^, is given by 


y K = E(x-y 1 ) K = 


(x-y^) 1 ^ f(x)dx. 


(7) 


These moments about the mean can sometimes be simplified if it 
is noted that 

E(Ex i ) = E[E(x i )3, (8) 

and 

E(cx) = c*E(x) . 


(9) 
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The second moment about the mean, a measure of the spread or 
dispersion of the distribution, is called the variance (y 2 ) and is 
obtained as follows 

y 2 = E(x-y 1 ) 2 (10) 

or ‘ 

Vg = E(x 2 ) 2y 1 E(x) + (y.^ 2 . 

But ... . ... .. 

(ll) 


( 12 ) 


(10a) 

The third moment about the mean, called the skewness (y ), is 

o 

a measure of the symmetry of the distribution and is given by 


---E(x) = 

and --- — : ... 

V -r 1 . E(x 2 ) = y 2 «. 

Therefore, . 

y 2 = y 2 '.- (y 1 ) 2 . 
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U 3 = E(x-y 1 ) 3 , <UL3) 

V *3 = ^ 3 * ~ 3 ^ 2 ,y l + 2y l 3, (13a) 

The coefficient of skewness, a^, is defined as 

a 3 “ Vi 3 /(P 2 ) 3 ^ 2 , (14) 

and is a measure of the skewness relative to the spread of the dis- 
tribution. If a > 0, the distribution is skewed to the right, 
o 

This means that a tail extends to the left. If < 0* the distri- 
bution is skewed to the left, and if = 0, the distribution is 
symmetrical . 

The fourth moment about the mean, called the kurtosis (y^), is 
a measure of the peakedness of the distribution, and is given by 

y 4 = E(x-y 1 > 4 , (15) 

V 4 = W 4 ' ~ 4 I j 3 '4 1 + 6y 2 '(y 1 ) 2 - 3(y 1 ) 4 . (15a) 

The coefficient of kurtosis is defined as 

= (h 4 )/(y 2 ^ 2, (16) 
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The value of a Ls to measure the peakedness of distributions 

relative to each other. The distribution with the largest is the 
most peaked. ' 
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cesnri b-jv.- - DISTRIBUTIONS APPLICABLE TO FATIGUE 

* 

u The statistical distributions which are frequently used in, 
and those which have been successfully applied -to, the explanation 
of fatigue phenomena are the following: 

' ~ - lV' - Gaussian or normal 


2. 

Lognormal 

3. 

Gamma 

4.- 

Erlangian 

5." 

Exponential 

6. 

: Extreme-value 

7. 

Weibull 

8. 

Beta 



The choice of the appropriate distribution to be applied to 
specific fatigue phenomena has to be based on the statistical, as 
well as, the phenomenological aspects of the generated fatigue data. 
These aspects are discussed next. 




15 


A. The Normal Distribution 

The normal is the most widely used of all distributions. It is 
best described as the distribution of events which are the result of 
the sum of many small effects . Many random events are not normally 
distributed, and therefore there would be no reason to expect a nor- 
mal distribution. However, one justification for examining the 
normal is the central limit theorem. 

The central limit theorem states that the distribution of the mean 
of n independent observations from any distribution with finite mean and 
variance approaches a normal distribution as the number of observations, 
n, approaches ro . This is an important principle, because 
although it applies to a large number, of observations , even a rela- 
tively small number of observations will tend to normality if the 
parent distribution does not deviate too far from the normal. 

The probability density function of the normal distribution is 
given by 

1 

f(x: M, °> = 



The mean is given by 


E(x) = 


g/2tt 


x exp 


h(^] 


2i 


I dx. 


( 18 ) 


or 


E(x) = p. 


( 18 a) 
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The variance of normal distribution is 

' ■ I ... - 

Var(x) = a^. (19) 

The two parameters of the normal distribution are the mean, y, and 
the standard deviation, o. The normal is symmetrical about the 
mean and is def ined for values of x between positive infinity and 
negative infinity. 

The cumulative density function of the normal distribution is 
given by 


FU) = 


This integral can only be approximately evaluated. Tables of its 
values exist for y = o and a = 1. 

For normal distribution = 0, which shows that the distribu- 
tion is symmetrical, an d ct 4 = 3 '(7, pp. 123-124). 
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B. The Lognormal Distribution 

A distribution closely associated with the normal is the lognor- 
mal. This is the distribution of a random variable whose logarithm 
follows the normal distribution. The lognormal probability density 
function is 


f(x: y, o) = 


xo 


"/2tt 


exp 




(log x-y")‘ 


( 21 ) 


x is defined from zero to positive infinity. The distribution is 
skewed to the right with skewness increasing as 0 increases. ' 

The mean of the lognormal is 


V = e , 


( 22 ) 


where 


n 


M" *■- l log X , 
1=1 


(22a) 


and 


.V log x i - v 

i=l 


(22b) 


and the variance is 



18 


2 2 

( e ° - 1 ). 


The coefficients of skewness 5 ' ^ :ur> ^ os ^ s are PP* 127-128) 


a 3 = 'if* (e°" + 2), 


and 


'2o" i 


2c"' 


= 3 + (e° - 1-y^ +. 3e + 6e 


+ 6 ), 



o 


(23) 


(24) 
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C. The Gamma Distribution 

The gamma distribution is useful for representing the distribu- 
tion of quantities which cannot be negative. It is appropriate to 
define the distribution of the times required for a total of exactly 
3 independent events to take place if they occur at a constant rate 
n. It could be used to represent times-to-failure of components if 
the subcomponents fail independently with a constant rate q. 

The gamma probability density function is 

< g 

f(x: 3, n) = ^ e“ nX v (25) 

It is defined for x 1 0, 3 > 0, and n > 0. r ( 3 ) is the gamma func- 
tion. 

The cumulative distribution of the gamma is 


fx 


F(x: 3, n) = 


_ n 


r(3) 


3-1 -nx , 
x e dx. 


(26) 


The 

accounts 

The 


gamma distribution has a wide variety of shapes , and this 
for much of the use of the model, 
mean is 




( 27 ) 


\ 


and the variance is 
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~ 'The- coefficients of skewness and kurtosis are (7, pp. 123-120). 


_2_ 

“3 ' ft * 


and 



3(6 + 2) 


3 
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D. The Erlangian Distribution 

! 

If n is restricted to positive integers, the gamma distribution 

\ 

is referred to as the Erlangian distribution. The distribution is 
sometimes more realistic in this form, as for most applications the 
fraction of an event has no meaning (7, p. 90). 
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E. .. -The Exponential. Distribution . • 

I 

The exponential distribution is often used as a times-to-failure 
distribution. It is used when the failure rate is assumed to be con- 
stant. It is the times-rto-failure distribution if. these failures are 
independent and happen at a constant average rate. 

The. probability . density function for the exponential distribu- 
tion is. (.7,. pp. 123-124)... 

. f(x:n) = ne _T1x , (29) 

and is defined for x > 0 and n > 0. The mean is 



( 32 ) 
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_ and' the: coefficient of kurtosis is 



p.tr r.v From the values of and it is seen that the exponential 
-distribution is skewed 'to the right and more peaked than the normal. 
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F. The Extreme-Value Distributions 

The extreme-value distributions should be considered in a dis- 
cussion of fatigue life distributions. In many applications the 
distribution of the largest or smallest elements of a sample are of 
interest. In failure analysis the distribution of the weakest com- 
ponents (smallest values) would be of interest. A distribution of 
the minimum of n independent values from a parent distribution that 
is unbounded to the left and is of exponential decreasing type is an 
extreme-value distribution. This distribution is the Type I for 
minimum values . 

In reality, failure times cannot be negative. Therefore, the 
real life distributions of times-to-failure should be bounded by 
zero on the left. One such extreme-value distribution is the Wei- 
bull, referred to as the Type III extreme-value distribution. 



25 


G. The Weibull Distribution 


The only distribution - that was actually devised for use with 
fatigue data is the Weibull. It is an extreme value distribution 
of the smallest values. Using the Weibull to represent the break- 
ing strength of a material has also been justified by Freudenthal 
and Gumbel (9). f , 

The probability density function for the Weibull is (7, pp. 
131-132) 


f(x. : 3, n, y) = — 


x ~ Y 


3-1 


exp [- 


x ~ Y 


(34) 


cr 3 > 0 is the shape parameter. 

Y > 0 is the location parameter, 
t) > 0 is the scale parameter. 

The mean of the function is (11, pp. 2-15) 


y = Y + n T 


(H 


The variance is (7, p. 132) 

a 2 = n 2 [r 


f| + xl 

- r 2 

fui] 

J 


J 


3 . 


(35) 


(36) 


The coefficient of skewness is (7, p. 132) 

_ r(i + 3/3) - 3 rd + 2/3) r(i + 1/3) + 2 [r(i + 1/3)] 3 


a. 


(r(i + 2/3) - Cr(i + i/ 3)] 2 > 3/2 


(37) 
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and the coefficient of kurtosis is (7, p. 132) 


_ r(i+4/6)-4r(i+3/s)[r(i+i/B)]+6r(i+2/e)[r(i+i/B)] 2 -3[r(i+i/B)] 4 . (38) 

' {r(i+2/B) - Cr(i+1/B)] 2 } 2 

The cumulative distribution function is given by 


F(x) = 1 - exp [- 



(39) 


This is the unreliability function. The reliability function is defined 
as 


R(x) = 1 - F(x), 


or 


R(x) = exp [- 


x-Y 
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l h J 




(40) 

(40a) 


and ^ 

RKT =ex P C 

If the location parameter is zero. 


x - y 


]. 


R(x) 


= exp 


x 

A 


(40b) 


then 


log 


r y 

1 


t 

X 

[R(x)J 


.4. 


and 


log log 


(rHoJ 


= $ log x - B log n* 


(40c) 
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This equation is of the form 

Y = AX + B 

Y “ log log (kc4t , 

X = log X, 

B = -g log n (41) 

Using Weibull probability paper and plotting the times-to- 
failure results, it can be determined if the Weibull distribution 
describes the data* If a straight line fits the plotted points, 
then the Weibull distribution is sufficient for describing the data* 

In addition, the parameters of the distribution can be determined 
from this plot. The slopeTof the line is B, the shape parameter, 
n is the abscissa corresponding to an ordinate value of 63.2 per- 
cent on the Weibull plot. If a straight line fits the plotted 
data points, then y = 0. If not, and a curve drawn through the 
points exhibits a concave downward behavior, then the data may be 
adjusted so that a straight line fits them satisfactorily. Then 
the value of y can be determined, by the method given by Lochner 
(18)/ :r ‘ / ‘ ~ " . 
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H. The Beta Distribution 

The beta distribution is useful to describe variations over a 
finite range. It has limited use for predicting cycles-to-failure 
but is included here for completeness. The beta probability density 
function is (7, pp. 127-128) 


f(x* 8 n) = ^ ^ * — 1 x n-1 (l _ 

n; r(8) r(n) u ' 


(42) 


for x defined over the interval zero to one, and n and y > 0. 
The mean is 


y = 


n 


8 + n 


(43) 


The variance is 


v2 = 


8 n 


(8 + n) 2 (8 + n + l) 

The coefficient of skewness is • 


(44) 


a = 2(8 - n) (n + 8 + 1) 
3 (8n) 1/2 (n + 8 + 2) 


1/2 


(45) 


The coefficient of kurtosis is 


_ 3(3 + n + 1) [2(ri + 3) 2 + 3n(3 + n - 6)] 
4 ~ 

3n(£ + n + 2) (3 + n + 3) 


( 46 ) 
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V THEORIES OF FATIGUE FAILURE AND THE DISTRIBUTIONS 
ASSOCIATED WITH THESE THEORIES 

In this section the three predominant theories of fatigue failure 
are discussed. Also an attempt is made to show how each one of these phe- 
nomenological theories leads to the use of a particular distribution. 

Many tests have been performed on fatigue failure, and some 
characteristics have been observed. Beyond a certain number of cycles 
of operation the occurrence of fracture is probabilistic (10, p. 214). 
This number is stress dependent. Many cases show a positive skewness 
to the frequency distribution of fatigue life (10, p. 214). The scat- 
ter in fatigue life decreases as the stress is increased (10, pp. 214- 
215)/ 

Fatigue life of a specimen can be divided into four (10, p. 218) 
stages. The first stage of n^ cycles is -the completion of work harden- 
ing. The second stage of n^ cycles is the time in which the first 
microcracks are formed. The third stage of n^ cycles is the period 
during which the submicrocracks grow and link to form a crack of detec- 
table size. During the fourth stage of n^ cycles these cracks propa- 
gate across grains until fracture or rupture occurs. The end of stage 
3 and the beginning of stage 4 are not clearly defined, and the period 
between formation of a detectable crack and rupture is a small part of 
the total life. Therefore, the contributors to the scatter in fatigue 
life are the second and third stages. 

It is the rate of crack growth which determines the number of 
cycles before failure after the first cracks have been formed. At low 
stress levels, just above endurance, cracks have been found to exist 
after 50% of fatigue life. At higher stress levels cracks appear just 
before failure (10, Figure 139, p. 208). 
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A-.,.: The . Weakest Link Theory ; 

The weakest link theory treats each component as a series of 
many subcomponents. It interprets the strength of the component in 
terms of the minimum values of the strengths of the subcomponents. 
Each link will have cracks with a certain distribution, and the 
component will fail when the weakest link fails. 

z~zz .:. The cracks or defects are distributed throughout the specimen. 

Of the n total subcomponents the least strong determines the 
strength. The distribution of interest will be that of the minimum 
values of the subcomponent strengths in n subcomponents. If the 
life cycles of the subcomponents have a probability density function, 
f (x) , the cumulative distribution function, F(x), is the proba- 

bility that the life cycles do not exceed x and is defined as 


F(x) 


[ f(x) dx. 


— 00 / 


(47) 


The life lengths of the aggregated component would be distributed 
according to the smallest order statistic-, thus. 


F n (x) = 1 - [1 - F(x)3 n . 


(48) 


Consequently, the probability density function, f n (x), of the 
smallest value of the n life cycles is (11, p. 2-6) 

f n (x) = n[l - F(x)] n_1 f(x). 


(48a) 



31 


~:The Weibull distribution is a smallest value distribution. If 
the subcomponents have a life, distribution of the Weibull form, then 


• 1 — F(x) = 1 - e . , (49) 

and the cumulative distribution of the aggregate from Equation 48 
becomes - . .acre.,.. — • 


..... ..... . r r _ 

F (xj = 1 exp -n — - _. (50) 

A . I h J J 

and the pdf is then given by (11, pp. 2-6) 




(51) 


There are also other subcomponent populations which will lead to 
the Weibull distribution for the specimen if the weakest link theory 
is accepted. Fisher and Tippet (12) have shown that many distribu- 
tions, including the normal, can have their smallest values 
distributed as an extreme “value distribution. 
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Br The Parallel Strand Theory 

To understand this theory it is convenient to consider the 
strands of a multistranded rope. The component is made up of subcom- 
ponents such as the strands of a rope and cannot fail until every 
strand has failed. Life lengths of the subcomponents determine the 
life pattern of the component. The life pattern of the component is 
the convolution of the life patterns of the subcomponents. 

If- the subcomponents had life lengths that were independent and . 
identically distributed, the life length for the component would be 
the n-fold convolution of the subcomponent distributions. If the 
subcomponents have a probability density function f(x), the life 
length pdf, g(x), of the component will be (11, pp. 2-7)' 

g(x) = [f ( x ) ] n< (52) 

where . 

[f(x)3 n * = [f(x)] *Cf(x)]^-D* (53) 

and 

Cf(x)]( n - 1 )*=[f(x)] *[f(x)] (n - 2 ^ (53a) 

and 

[f(x)3 2 * = [f(x)] *[f(x)] ' (54) 


b 
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or 


Cf(x)3 


2 * 


f(t) f(x-t) dt 


(54a) 


Assume the life distributions of the parallel strands are expo- 
nentially distributed; then 


and 


and 


or 


f (x) = ne 


-nx 


(55) 


[f(x)3 


2 * 


r x 




-T)t 

-n(x-t) 


ne 

ne 


k * 

. W 4 


dt 


(56) 


Cf(x)] 


2* 2 ? 
= n 

0 J 


e“ nx dt 


(56a) 


r-j-, -.-,2* 2 -nx 

Lf (x) J = n t e 


(56b) 


Therefore , 
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r - f . -,2* i 2 -rpi 

[f(x)3 = n x e 


(56c) 


Similarly , 


rx . t \ -.3* 3 

[f(x)3 = n 


rx 


t e _Ilt e-" (x - t) dt. 


(57) 


[f(x)3 3 ** = x 2 e~ nx , 


(57a) 


and (11, pp. 2-7 and 2-8) 


rj=r v-,4* 14 3 -nx 

Cf(x)3 = gj n x e 


(58) 


Consequently, 


[f(x)3 n> = ( R _3j! n n x irl e -r|X j 


(59) 


or 


g (x) = [f(x)] 
n 


n n-1 
_ n x -nx 


r(n) 


(59a) 


It is seen from Equation 59a that the pdf of the component is defined 
by a gamma distribution. 

If the number of strands approaches infinity, as it would for a 
metal specimen, the life length distribution g^(x) approaches the 
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normal function. As the shape parameter of the distribution, n, 

• i 

increases, the gamma distribution tends to normality with a mean of 

2 

n/ri and a variance of n/n (11, p. 2-8). 
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C. The Proportional Effect Theory 

i - • 

A component is assumed to fail when the size of a crack reduces 

the cross sectional area to a certain value. A crack propagates at 

an exponential rate. The crack length is proportional to the length 

of the preceeding stage. Assume the size of a fatigue crack at various 

stages of its growth can be represented by the sequence x-j_ < < ... 

x < . . . < x , where x is the size of the crack at the r^ a stage . 
r n’ r 6 

x is the size of the crack when the cross sectional area is reduced 
n 

to a value that cannot sustain the applied stress , and rupture occurs . 

"th 

The crack growth x^ - x^_^ at the i stage is proportional to the 
crack size x^_^ of the preceeding stage, or 



9 .x. . 
i l-l 


(60) 


x^ can be interpreted to be the size of minute flaws in the original 
component. 9^, 9^, ... are independently distributed proportionality 
constants (11, p. 2-8). Then 


x. 

l 


9.x. , + x . . , 
i l-l l-l’ 


and 


(60a) 


x i = x i-i (a i + 1) ' 


(60b) 


x i-l = x i-2 (8 i-l + 1) * 


(60c) 


It follows that 
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X 2 = X l (3 2 + 1 ^» 


(60d) 


-c-_.no::... ; y :ix-' =--x (3 + 1), 

l o 1 




-Then 


(60e) 


^ ^ ^ X o ( V 1)(8 2 +1)5 


(61) 


" = " X q( ^ 1 + 1) ( ••• (61a) 


and 


x; = x (3 n +l)(3 0 +l) ... (3. n +l)(3.+l). (61b) 

l . . o 1 2 l-l l 


■c:j“:*:r,The component is assumed to fail when the size , of the crack 

reaches x^. The characteristic life length of the component is the 

distribution of x , where 
n 




x = x (3_ +1)(3 0 +1) .... (3 +1). (61c) 

n o 1 2 n 


From Equation 61c it may be seen that x n is the product of indepen- 
dently distributed random variables. The logarithm of x^ is the sum 
of independent random variables , or 


log x = log x + log (3+1) + ... + log (3 +1). 
n o 1 . n 


( 62 ) 
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The central-limit theorem states that log x^ is approximately nor- 
mally distributed if n is iarge. If log x n is normally distributed, 
then x^ is lognormally distributed (11, p. 2-8). 

At higher amplitudes of alternating stress many cracks are formed. 
The growth rate of each of these cracks follows the proportional-effect 
theory. However, the crack lengths interlink and cause a random reduc- 
tion of component area. At lower stress levels only one crack usually 
exists. The growth of this single crack also follows the proportional- 
effect theory. 

. Freudenthal (13) proposes a derivation for the distribution of 
component failure using the proportional-effect theory. It is based on 
the assumption that a single crack is formed and it propagates to cause 
failure. 

Let c , c , . . . , c, be consecutive cycles applied to a specimen at 
a constant stress amplitude, and let the extent of damage done to the 
area M by the cycles c^, . • . , c^ be M^. The increase in the damage 

done by cycle c^ is This increase is proportional to the 

effect of the cycle c^_^ and is related by some function 4>(M^); conse- 
quently. 


and 


- “k-i '= C k+ ( V’ 


(63) 


C, = 


M, - M, , 
k k-1 


k <j>(M k ) 


(64) 


Ji ° k = 


(65) 
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area , 


If each cycle contributes only slightly to the 
the sura can be replaced by 



disruption of the 



r** 


n 


M > 
o 


dM 

<i>(M) * 


( 66 ) 


If the effect of each cycle is directly proportional to the extent of 
damage produced by the previous cycle, then ( M ) is constant, and 


I 


k=l 


(M 

n 


dM _ 
kM " 


M * 
o 


log M /M 
• n o 


(67) 


where M is the extent of the initial damage, 
o 

The central limit theorem states that the sum of n independent 
ramdom variables tends toward a normal distribution as n increases. 
If the number of cycles is large , then log M^/M^ is normally distri- 


buted. 

Making the assumption that the average rate of damage is pro- 
portional to the damage produced by a given number of cycles and 
that the cycles to produce a given amount of damage is inversely 
proportional to the rate, the number of cycles to fracture is 
inversely proportional to M. The reciprocal transformation of the 
lognormal distribution is also lognormal, log x = - log 1/x. The 
distribution of fatigue life at a given stress level is, therefore, 
lognormal (10, pp. 235-236). 

The derivation of Freudenthal, which shows the life length dis- 
tribution of the component to be lognormal, explains the previously 
mentioned positive skewness of failure distributions. The interlinking 
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of cracks to cause random reduction in the cross sectional area at 
high stress levels indicates that these life length distributions are 
normal. This has been supported by our results and the results of 
other testing programs. 
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VI” "ANALYSIS OF THE FATIGUE RELIABILITY AND STATIC STRENGTH DATA 

GENERATED. IN THIS RESEARCH 

The decision as- to -'which distribution should be chosen must be 
based on the plots of the cycles-to-failure and stress-to- 
f allure data, and phenomenological reasoning. The data was generated 
on the three combined bending- torsion fatigue reliability research 
machines designed and built at The University of Arizona. The times- 
to-failure were determined to the nearest second by precision clocks* 
which started' as soon as -a -test was underway and stopped when a micro- 
limit switch- cut the power to the clock off as soon as the specimen 
failed. The r time was then converted to cycles-to-failure using the 
rotational speed of the machines. The times-to-failure and their 
conversion to cycles-to-failure are given in Appendix A. The speed 
of” the machines was calibrated to +_ 5 rpm. The rotational speeds of 
the machines are:" ' tv; - '- 

-iz-.-ac.-v- Machine 1 - - - 1786 +_ 5 rpm 

Machine 2 ■- 1784 _+ 5 rpm - 

Machine 3 . 1780 + 5 rpm 


For the data generated in this research and given in Tables 1, 
2,” and 3, it was decided to first plot the cycles-to-failure and 
stress-to-failure data on probability paper, study the results, cor- 
relate them -with the previously discussed behaviors of the statis- 
tical distributions and the phenomenological aspects , and then draw 
conclusions as to which distribution best represents the data. 

Three stress levels of stress ratio 00 were run with a positively 
driven revolution counter. 



42 


Based on the discussion of the distributions relevant to fatigue 
and static strength, it was decided to use only the normal, lognormal 
and Weibull distributions. Consequently, all of the cycles-to-failure 
data were plotted against median ranks ( 17, Table 1 ) on Weibull 

and lognormal probability paper for all stress lev- 
els of stress ratios ~ and 0.70. The data for the two highest stress 
levels at each stress ratio were plotted against median ranks on 
normal probability paper also. • 

In Figure 5 are plotted the cycles-to-failure data for r = » 

s 

on Weibull probability paper. In Figure 6 are plotted the cycles-to- 

failure data for r = 00 and s = 154,000 psi and 121,500 psi, on nor- 

mal probability paper. In Figure 7 are plotted the cycles-to-failure 

data for r = °° and s = 154,000 psi, 121,500 psi, and 104,500 psi 

on lognormal probability paper. In Figure 8 are plotted the cycles- 

to-failure data for r = » and s = 86,000 psi and 78,000 psi on log- 

s a 

normal probability paper. 

In Figure 9 are plotted the cycles-to-failure data for r = 0.70 

s 

on Weibull probability paper. In Figure 10 are plotted the cycles- 
to-failure data for r =0.70 and s = 110,500 psi and 97,500 psi on 
normal probability paper-. . In Figure 11 are plotted the cycles-to- 
failure data for r = 0.70 and s = 110,500 psi, 97,500 psi, 76,500 
psi, and 70,000 psi on lognormal probability paper. Lastly, in Figure 
12, are plotted the static tensile ultimate strength data for the 
notched specimens on normal probability paper. 

In Figure 5, the data for the stress level of 154,000 psi plots 
concave upward. The location parameter, y, could not be adjusted to 
give a better straight line fit ( 18, pp. 5-8 ). The raw data would 
lend themselves to a fairly good straight line fit at the stress levels 
of 121,500 psi, 104,500 psi, 86,000 psi, and 78,000 psi, if these lines 
were drawn. The lines on this Weibull probability paper were not 
drawn so that the raw data would stand out. 
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Iri r Figure 6, the data for the stress level of 154,000 psi is 
concave upward' indicating a distribution that is skewed to the left 
(16 i p. 14' ). " The data for the'stress level of 121,500 psi gives 
a Straight line fit with good correlation, as shown in Table 4, on 
normal paper. •' ' 

Figures 7 and 8, and Table 4 show that the data for the stress 
levels of 121,500 psi, 104,500 psi, 86,000 psi, and 70,000 psi give 
straight line fits with good correlation on lognormal paper . A 
straight line can be fitted to the data for the stress level of 
154,000 psi, although with not as good a correlation as at the 
other stress levels /The lack of an extremely good fit may be ex- 
plained by ‘the fact that only twelve data points are available at 
this stress level. ‘ : 1 ' 

* Figure 9, where the plots of the data for stress ratio of 0.70 
are given on W.eibull probability paper, shows that a straight line 
would fit the data at ^he lower stress levels of 76,000 psi and 70,000 
psi with good correlation. The data for the higher stress levels of 
110,500 psi and 97,500 psi, are concave downward. The adjustment of 
the location parameter did not yield a better straight line. The 
attempt to adjust for better fit was hampered by the lack of a larger 
number of data points. . 

Figure 10, where the data for the stress ratio of 0.70 are plot- 
ted on normal probability paper, indicates concave downward curves 
for the stress levels of 110,500 psi and 97,500 psi. 

Figure 11 and Table 4 show good straight line fits at all stress 
levels for a stress ratio of 0.70 on lognormal probability paper. 

In Figure 12 and Table 4, the ultimate strength data exhibit a good 
straight line characteristic on normal probability paper. The line 
has a steep slope indicating a narrow spread of the data. This 
is to be expected. Static strength distributions 

are usually normal. Juvinall ( 6, p. 351 ) states that static tests 
have a small statistical variation. Bompas-Smith ( 14, p. 344 ) also 
states that the probability density function of tensile tests can be 
expected to be normal. 



Having described each figure, an overall analysis of each will 
be made in conjunction with the phenomenological aspects of fatigue, 
discussed previously. A straight line can be fitted to all the log- 
normal plots with fairly good correlation, as may be seen in Figures 
7, 8, and 11, and Table 4. Phenomenologically, the acceptance of 
the proportional- effect theory, discussed earlier in this report, 
would result in the acceptance of the lognormal as the life length 
(cycles-to-failure) distribution of components subjected to fatigue. 

The conclusion that the characteristics of the lognormal distribution 
are associated with fatigue failures and that the lognormal is the 
failure governing distribution for fatigue is supported by b oth theory 
and experimental results. 

Herd ( 19, p. 5 ) also reasons that the lognormal is an appro- 
priate distribution for the cycles-to-failure data. He states that the 
lognormal distribution applies to situations in which several indepen- 
dent factors influence the outcome of an event, not additively, but 
according to the magnitude of the factor and the age of the item at 
the time the factor is applied. If the effect of each impulse is 
directly proportional to the momentary age, x, of the item, then log x 
would be normally distributed. Consequently, the x's would be log- 
normally distributed. 

Yokobori ( 9, p. 194 ) states that a positive skewness to the 
distribution of fatigue life often exists , and the logarithms of 
cycles-to-failure can be approximated by a normal distribution. 

Results of tests that show the positive skewness are given by Yoko- 
bori ( 10, pp. 211 - 212 ). The derivation by Freudenthal ( 13 ) 
and the discussion of the proportional-effect theory ( 11, pp. 2-8 ) 
also result in the lognormal as the distribution of cycles-to-failure 
for fatigue. Bompas-Smith ( 14, p. 345 ) states that fatigue results 
at a constant stress level frequently conform to a lognormal distri- 
bution. F. Epremian and R. F. Mehl ( 15 ) suggest the values of the 
logarithms of cycles-to-failure are normally distributed about a 
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mean value. Juvinall ( 6, pp. 350-351 ) shows results of tests of 
fatigue life data that approximate the lognormal distribution. 

Figure 10 shows curves that are concave downward. James R. King 
( 16, p. 7 ) states that a concave plot on normal probability paper 
indicates a right-skewed distribution and that a logical choice would 
be the lognormal. Bompas-Smith ( 14, Figure 12, p. 349, and Figure 
15, p. 350 ) confirms that a curve of this shape on normal probability 
paper gives rise to a straight line fit on lognormal probability 
paper. 

Figure 9 indicates that the lognormal distribution would provide 
good fits at stress levels of 110,500 psi and 97,500 psi for the stress 
ratio of 0.70 because the curves are concave downward and Bompas-Smith 
( 14, Figure 12, p. 349 and Figure 15, p. 350 ) shows that plots of 
this shape on Weibull probability paper give a straight line on log- 
normal probability paper. 

The extreme-value function could phenomenologically be the life 
length distribution if only the weakest of the specimens were tested 
to failure. This does not completely describe the testing program 
by which this data was generated. The specimens to be tested were 
randomly selected and all were tested to failure. Although the 
Weibull plots show that a straight line fit to much of the data 
is plausible, it is in no case better than the straight line fit pro- 
vided by the lognormal plots . 

It is difficult with the eighteen data points at each stress 
level to determine the exact shape of the Weibull plots . It would 
be impossible to determine the shape with the five or seven data 
points recommended by King ( 16, p. 12 ). This is true for any 
other distribution. With only a few data points, no absolute state- 
ments can be made from the plots. That is why the plots must be 
used along with phenomenological reasoning to determine the failure 
governing distributions. 
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Phenomenological reasoning can also be used to justify the use 
of the normal as the failure distribution at the high stress, or 
low fatigue life, levels. Juvinall ( 6, p. 351 ) states that short 
life fatigue tests approach static test characteristics which have 
small statistical variations. Bompas-Smith ( 14, p. 344 ) states 
that if the strength of a component is a function of several varia- 
bles, the failure distribution tends to normality. Yokobori ( 10, 
p. 215 ) states that the scatter of fatigue life increases as the 
stress level decreases. This is verified by the data in Tables 1 
and 4, and Figures 6 and 10. 

A computer program was written to fit the best straight line to 
the data ( straight cycles-to-failure, and logarithms of the cycles 
to-failure). This program uses the least squares method and fits the 
best straight line to the data for stress ratios <=° and 0.70 on normal 
and lognormal plots. The program also fits the best straight line 
to the data for the stress ratio of 0 on the normal plot and computes 
the correlation coefficient. 

The correlation coefficients for the data fits at the various 
stress levels differ only slightly between the normal and lognormal 
distributions. The maximum difference is 2.2%, as may be seen from 
Table 4. Basing the analysis on the relatively few data points makes 
it impossible to discriminate between the normal and the lognormal 
distributions solely on the basis of the correlation coefficient. 

The straight line fit to the data for stress ratio 0 is very 
good, because a very high correlation coefficient was obtained. The 
coefficients of correlation and the equations of the best straight 
lines are given in Table 4. Phenomenological reasoning, experimental 
results and graphical analysis dictate that the normal distribution 
should best describe the stress -to-failure data a stress ratio of 0. 

Although in many cases the straight line fit to the data on 
normal probability paper gives a higher correlation coefficient, 
the lognormal has been chosen to represent the cycles-to-failure 
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distribution. The small sample size used in this research has not 
provided sufficient opportunity for significant discrimination be- 
tween the normal and the lognormal distributions. Furthermore, the 
lognormal distribution has been phenomenologically justified and is 
seen to better represent the data, when all stress levels are consid- 
ered, than. any other applicable distribution discussed in this re- 
port. . . * 
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VII CONCLUSIONS 

1. The Weibull distribution fit to the cycles-to-failure data 
at various stress levels ( Figures 5 and 9 ')' shows that the distri- 
bution might approximate the cycles-to-failure distribution of the 
specimens; however, the data points do not appear to lend themselves 
to a good straight line fit because of the kinks in the plot. Sta- 
tistically, the Weibull distribution is an extreme-value distribution 
and does not describe the type of data generated in this research 
program. 

2. The normal distribution fit to the cycles-to-failure data 
at the various stress levels for the stress ratio of 00 and 0.70 

( Figures 6 and 10, and Table 4 ) shows that the normal distribution 
might represent the cycles-to-failure distributions of the specimens , 
because there is good straight line fit and there is no significant 
difference between the correlation coefficients for the normal and 
the next appropriate distribution, the lognormal. Phenomenologically, 
the normal distribution can be justified to approximate the cycles- 
to-failure distribution at the highest stress levels only . 

3. The lognormal distribution fit to the cycles-to-failure 
data ( Figures 7, 8, and 11, and Table 4 ) is good at all stress 
levels for stress ratios of » and 0.70. The lognormal distribution 
phenomenologically describes the cycles-to-failure distributions of 
the specimens best at all stress levels . 

4. Phenomenologically, statistically and probabilistic-graph- 
ically the normal distribution gives the b est fit to the static ul- 
timate strength data for the stress ratio of 0 ( Figure 12 and 
Table 4 ). The straight line fits the data on normal probability 
paper with good correlation and has a steep slope showing a small 
dispersion. 
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5. The difficulty of discriminating between the normal and log- 
normal distribution fits to the cycles-to-failure data is attributed 
to the small sample size; namely, 12 specimens at the highest stress 
levels and 18 specimens at the other stress levels. It would be of 
great interest to see the degree of discrimination achieved when 
sample sizes of 50 or more are tested at each stress level. 

6. The fact that the correlation coefficients for the straight 
line fits to the static ultimate strength data for notched specimens 
on the normal and the lognormal bases have no significant difference 
( Table 4 ) may again be attributed to the small sample size tested, 
namely 10. 

7. The phenomenological reasoning leads to the conclusion 
that, were the cycles-to-failure data those of field failures or only 
of the failures from a larger sample, all of which were not tested to 
failure, the Weibull distribution would be the most appropriate dis- 
tribution to represent such data. The primary reason for this is 
that such data would be the failures of the weakest of such components 
in field operation or in test, consequently , ^conf o rming to the extreme- 
va lue distrib ution theory. The cycles-to-failure data generated in 
this research are those of the whole sample being tested to failure; 
therefore, the data is that of the weakest, as well as, of the strong- 
est specimens failing, hence not conforming to the extreme-value theory 
represented by the Weibull distribution. 

8. Phenomenological reasoning leads to the conclusion that the 

cycles-to-failure data at all stress levels would be best represented 
by the lognormal distribution, and,' except at the highest alternating 
bending stress levels, the lognormal distribution should be used ex- 
clusively for the cycles-to-failure data of the type generated in this 
research. - 

9. Phenomenological reasoning also leads to the conclusion that 
at the highest alternating bending stress levels, the normal distri- 
bution can be used to approximate the lognormal distribution. This 
provides a computational advantage when calculating the reliability 
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of a component by the design-by-reliability methodology. 

10. A conclusion of caution is in order when attempts are made 
to apply the previous conclusions to test conditions not identical, 
or closely related, to those used in the generation of the data for 
this research. More complex test and field loadings will cause fail- 
ures not represented. by any one of the following idealized theories: 
he weakest link theory, the parallel strand theory and the propor- 
tional effect theory. Under these conditions some combination of 
these theories would be in effect. The cycles or times-to-failure 
data would then exhibit complex behavior not representable by any one 
distribution discussed in this report. 
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VIII. RECOMMENDATIONS 

1 v There is not enough statistical evidence to discriminate 
between the normal and the lognormal at most stress levels. This 
is the result of testing 12 or 18 specimens at each stress level, 
and reflects the need to test more specimens at each stress level: 
preferably 50 or more . 

2. A test program should be initiated to test the weakest com- 
ponents of a population. and„ examine the cycles-to-failure behavior 
using the extreme-value distributions. These life length distribu- 
tions would be of interest to design engineers because it is the 
weakest parts that fail in actual service. 

3. Research of the type leading to this and the previous three 
reports should b e continued to acquire the vast data needed for the 
effective application of the design-by-reliability methodology. 

4. Statistical distributions which may represent the more 
complex test, or field 3 loading situations should be developed and 
studied. 
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FIGURE L CYCLES-TO- FAILURE DISTRIBUTIONS AND ENDURANCE STRENGTH DISTRIBUTION FOR STRESS RATIO OF 



160 

150 



FIGURE 2. CYCLES-TO-FAILURE DISTRIBUTIONS FOR A STRESS RATIO OF 0.70 AND ENDURANCE STRENGTH DISTRIBUTION 
FOR A STRESS RATIO OF 0.90. 


Stress- to- failure 


57 



Cycles-to- failure, N 


FIGURE 3. DISTRIBUTIONS OF THE CYCLES-TO-FAILURE AND STRESS-TO- 

FAILURE AT A SPECIFIC LIFE (10 5 CYCLES IN THIS EXAMPLE) 
AS FOUND BY FATIGUE TESTS AT A CONSTANT STRESS RATIO. ' 






59 






v 

u 

0 

rH 

•H 

fO 


W 

(D 

•H 

V 


o 

o 

O 

o 

o 

O 

O 

O 

o 

o 

O 

• 

• 

• 

• 

• 

* 

♦ 

• 

• 

• 

* 

CD 

in 

o 

o 

o 

o 

o 

in 

CO 

CN 

H 

cn 

CD 

av 


in 

CN 







aanxypj 


Xuaoaaj 


FIGURE 5. CYCLES-TO- FAILURE DATA FOR STRESS RATIO OF ® PLOTTED ON WEIBULL PROBABILITY PAPER. 
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FIGURE 7. CYCLES-TO-FAILURE DATA FOR BENDING STRESSES OF 154,000 PSI, 121,500 PSI 
104,500 PSI AT STRESS RATIO OF «> PLOTTED ON LOGNORMAL PROBABILITY PAPER 




FIGURE 8. CYCLES-TO-FAILURE DATA FOR BENDING STRESSES OF 86,000 PSI AND 78,000 PS I AT STRESS 
RATIO OF « PLOTTED ON LOGNORMAL PROBABILITY PAPER. 










ncuni: 3. CTCtXS-TO-rAILUKE DATA FOR STRESS RATIO OF 0.70 PLOTTED ON WEI BULL PROBABILITY PAPER 







95.0 



FIGURE 10. CYCLES-TO-FAILURE DATA FOR BENDING STRESSES OF 110,500 PSI AND 97,500 PSI AT A STRESS 
RATIO OF 0.70 PLOTTED ON NORMAL PROBABILITY PAPER. 
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Cy cles - t o- f a i lure 

FIGURE 11. CYCLES-TO-FAILURE DATA FOR BENDING STRESSES OF 110,500 PSI , 97,500 PSI, 76,500 PSI 
AND 70,000 PSI AT STRESS RATIO 0.70 PLOTTED ON LOGNORMAL PROBABILITY PAPER. 









FIGURE 12. STATIC TENSILE ULTIMATE STRENGTH DATA FOR NOTCHED SPECIMENS PLOTTED ON NORMAL 
PROBABILITY PAPER. 




CYCLES -TO -FAILURE FATIGUE DATA FOR SAE 4340 STEEL 
MIL-S-5000B, CONDITION C4 ROCKWELL C 35/40 
FOR STRESS RATIOS OF co and 0.70 
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*The variability is the standard deviation of the distribution about the mean. 

**The coefficient of variation is the ratio of the variability to the mean expressed as a percentage 
***A11 cycles- to- failure rounded to the nearest 50 cycles. 

**^*A11 stresses rounded to the nearest 500 psi. 
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TABLE 2 

STATIC ULTIMATE AND BREAKING STRENGTH DATA AND 
RESULTS FOR NOTCHED SPECIMENS* 

(Stress ratio «* 0) 


Test 

No. 

Ultimate Load 
1000 lbs. 

Breaking Load 
1000 lbs. 

Ultimate 
Strength 
psi. ** 

Breaking 
Strength 
psi. ** 

1 

49.3 

47.0 

253,500 

305,000 

2 

49.6 

47.0 

255,000 

305,000 

3 

49.4 

46.3 

254,000 

299,500 

4 

50.3 

47.4 

259,000 

299,500 

5 - 

48.8 

46.0 

251,000 

306,500 

• 6 

49.2 

46.0 

• 253,000 

302,500 

7 

49.6 

46.8 

255,000 

304,500 

8 

49.8 

47.1 

256,000 

305,500 

9 

50.5 

47.7 

260,000 

309,500 

10 

49.9 

47.5 

256,500 

302,000 


Normal Distribution Parameters 
of Ultimate Strength of Notched 
Specimens: 

Mean » Sy n » 255,500 psi 

Standard Deviation ■ 6L “2 

S Un 


Normal Distribution Para- 
meters of Breaking Strength 
of Notched Specimens: 

Mean « 304,000 psi 

500 psi Standard Deviation «* 

GZ 3,000 psi 


*Specimen diameter at the base of the notch is 0.4975 in. which 
gives an area of 0.1944 sq. in. 

**A11 strengths rounded to nearest 500 psi. 




STATIC YIELD, ULTIMATE AND BREAKING STRENGTH DATA AND 
RESULTS FOR UNNOTCHED SPECIMENS 
(Stress ratio =0) 
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APPENDIX A-l 


CONVERSION OF TIMES-TO-FAILURE DATA TO CYCLES-TO-FAILURE 
FOR 154,000 PSI ALTERNATING STRESS AND STRESS RATIO “ 


Item 

No. 

Spec. 

No. 

Machine 

No. 

Time to 
Failure 
hr :min:sec 

Cycles to 
Failure 

Log e 

Cycles 

Median 
Ranks * 

1 

341 

1 

0:00:56 

1£67 

7.4188 

5.613 

2 

365 

1 

0:01:05 

1,935 

7.5678 

13.598 

3 

342 

1 

0:01:23 

2,471 

7.8124 

21.669 

4 

193 

2 

0:01:33 

2265 

7,9248 

29.758 

5 

196 

2 

0:01:37 

2^84 

7.9669 

27.853 

6 

204 

2 

0:01:40 

2,973 

7.9973 

45.951 

. 7 

166 

2 

0:01:40 

2,973 

7.9973 

54.049 

8 

133 

2 

0:01:42 

3,033 

8.0173 

62.147 

9 

225 

2 

0:01:43 

3063 

8.0271 

70.242 

10 

220 

2 

0:01:47 

3182 

8.0653 

78.331 

11 

191 

2 

0:01:47 

3182 

8.0653 

86.402 

12 

163 

2 

. 0:01:50 

3271 

8.0929 

94.387 


"(17, Table I) 



APPENDIX A-2 


CONVERSION OF TIMES TO FAILURE DATA TO CYCLES-TO-FAILURE 
FOR 121,500 PSI ALTERNATING STRESS AND STRESS RATIO « 


Spec. 

No. 

Machine 

No. 

. Time to 
Failure 

Cycles to 
Failure 

h°g e 

Cycles 

Median 

Ranks 

96 

1 - 

i 

3 

o 

7,112 

8.8695 

3.778 

86 

1 

> 

<y 

U 

7,622 

8.9388 

9.151 

69 

1 

a 

a) 

> 

7,717 

8.9512 

14.581 

89 

1 

U 

8,015 

8.9891 

20.024 

22 

. 1 

H 

8,088 

8.9981 

25.471 

20 

1 

a) 

> • 

*H 

8,376 

9.0331 

30.921 

76 

1 

•P 

•H 

8,860 

9.0893 

36.371 

111 

1 

w 

o 

8,925 

9.0966 

41.828 

101 

1 

A 

■p 

%092 

9.1151 

47.274 

73 

1 

•H 

£ 

>> 

9,?61 

9.1336 

52.726 

64 

1 

H 

-P 

9,302 

9.1380 

58.177 

80 

1 

O 

<u 

£ 

9,362 

9.1444 

63.629 

52 

1 

nd 

9,747 

9.1874 

69.079 

59 

1 

Q) 

B 4 
3 u 

9818 

9.1920 

74.529 

109 

1 

W 0) 
d -p 
a> c 

9990 

9.2093 

79.976 

82 

1 

6 2 
o 

w a 

10347 

9.2444 

85.419 

42 

1 

0) 

H fi 

o o 

19353 

9.2450 

90.849 

63 

1 

*>-> *H 
o -p 

1Q540 

9.2630 

96.222 


APPENDIX A-3 


CONVERSION OF TIME-TO-FAILURE DATA TO CYCLES-TO-FAILURE 



FOR 104,500 

PSI ALTERNATING 

STRESS AND 

STRESS RATIO 

OO 

Spec. 

Machine 

Time to 

Cycles to 

Log 

Median 

No. 

No. 

Failure 

Failure 

Cycles 

Ranks 

39 

1 

1 

3 

H 

o 

16^58 

9.6963 

3.778 

55 

1 

> 

0J 

u 

16£00 

9.7111 

9.151 

99 

2 

p 

<u 

16^20 

9.7362 

14.581 

33 

1 

£ 

£ 

18^36 

9.8056 

20.024 

71 

1 

H 

19352 

9.8075 

25.471 

65 

2 

<u 

> 

2Q376 

9.9319 

30.921 

46 

2 

•H 

P 

•H 

21,080 

9.9561 

36.371 

85 

2 

CO 

O 

04 

213.92 

9.9614 

41.823 

113 

1 

X 

p 

21^204 

9.9619 

47.274 

26 

1 

•H 

£ 

>> 

22,544 

10.0232 

52.726 

49 

2 

rH 

P 

2^886 

10.0383 

58.177 

110 

1 

O 

<D 

A L‘ 1 

2^640 

10.0707 

63.629 

82 

2 

• l-J 

r O 

24304 

10.0984 

69.079 

13 

2 

a) 

25^96 

10.1344 

74.529 

58 

1 

to a) 
rtf p 
<D 0 

27,024 

10.2045 

79.976 

88 

1 

6 o 

CO O 

27250 

10.2128 

85.419 

62 

2 

a) 

rH P 

o o 

27:464 

10.2206 

90.849 

104 

2 

t>> *H 
O P 

27,558 

10.2241 

96.222 



APPENDIX A-4 


' CONVERSION OF TIMES-TO-FAILURE DATA TO CYCLES-TO-FAILURE 



FOR 86,000 

PSI ALTERNATING 

STRESS AND 

STRESS RATIO 

OO 

Spec. 

Machine 

Time to 

Cycles to 

L °g e 

Median 

No. 

No. 

Failure 

Failure 

Cycles 

Ranks 

45 

1 

o 

60,002 

11.0021 

3.778 

94 

1 

> 

<u 

u 

64,857 

11.0799 

9.151 

102 

1 

c 

<d 

64,997 

11.0821 

14.581 

118 

1 

> 

•H 

U 

67,757 

11.1237 

20.024 

114 

1 


71^.66 

11.1728 

25.471 

29 

1 

iH 

<D 

> 

71,556 

11.1782 

30.921 

35 

1 

•P 

•H 

71,968 

11.1840 

36.371 

18 

1 

W 

o 

P4 

7^089 

11.1994 

41.823 

25 

1 


74,690 

11.2211 

47.274 

30 

. 1 

•H 

£ 

74,698 

11.2212 

52.726 

31 

1 

i — 1 

+J 

75,662 

11.2340 

58.177 

34 

1 

O 

a) 

u 

78,852 

11.2753 

. 63.629 

70 

1 

•rH 

•x) 

X) * ' . 

82^48 

11.3175 

69.079 

75 

1 

<U u 

3 S 

w a 
rd P 
<D O 

83,812 

. 11.3363 

74.529 

66 

1 

8^349 

11.3662 

79.976 

48 

1 

6 o 
w £ 

96,004 

11.4722 

85.419 

19 

1 

a) o 

* — i *H 

98,468 

11.4975 

90.849 

81 

1 

O -H 
>. P 

O i— f 

107,415 

11.5845 

96.222 



APPENDIX A-5 

^CONVERSION OF TIMES-TO -FAILURE DATA TO CYCLES-TO-FAILURE 



FOR 78,000 

PSI ALTERNATING 

STRESS AND 

STRESS RATIO 

CO 

Spec. 

L MacKine 

I Time to 

Cycles to 

L°g e 

Median 

No.' 

No. 

Failure 

Failure 

Cycles 

Ranks 

215 

‘2 

0:52:18 

93,303 

’ 11.4436 

3.718 

167 

~2 

0:57:58 

103,413 

11.5465 

9.151 

121 

. *2 

1:09:35 

124,137 

11.7292 

14.581 

135 

~2 

'1:10:36 

125,950 

11.7436 

20.024 

145 

"2 

1:11:38 

127,794 

11.7582 

25.471 

198 

-2 

. 1:21:43 

145,783 

11.8899 

30.921 

128 

"2 

‘1:26:59 

155,178 

11.9523 

36.371 

156 

;2 

1:36:34 

172,275 

12.0569 

41.823 

171 

~2 

1:36:44 

172,572 

12.0586 

47.274 

123 

2 

1:36:47 

172,662 

12.0591 

52.726 

159 

2 

'1:39:30 

177,508 

12.0868 

58.177 

134 

' 2 

1:39:53 

178,192 

12.0906 

63.629 

~183 

2 

1742 : 30 

182,860 

12.1165 

69.079 

186 

2 

1:42:45 

183,306 

12.1189 

74.529 

214 

2 

1:49:46 

195,824 

12.1850 

79.976 

125 

2 

• 1:50:22 

196,894 

12.1904 

85.419 

148 

2 

1:53:48 

203,019 

12.2211 

90.849 

219 

2 

1:54:56 7 

205,041 

12.2310 

96.222 



I APPENDIX A-6 

CONVERSION OF TIMES -TO -FA I LURE DATA TO CYCLES-TO-FAILURE 
"FOR 110 ,500 PSI ALTERNATING STRESS AND STRESS RATIO 0.70 



' Spec. 
No. 

Machine 

--No. 

Time to 
^Failure 

Cycles to 
Failure 

L°g e 

Cycles 

Median 

Ranks 

=397 ” 


~"0T03 : 01 

5,388 

8.5919 

5.613 

' 395 

I 

0:03:09 

5,626 

8.6351 

13.598 

^ 432 

i 

0:03:18 . 

5,894 

8.6817 

21.669 

"387 

I 

0:03:23 

6,043 

8.7067 

29.758 

. ‘ 335 

i 

0:03:24 

6>072 

8.7114 

37.853 

' 358 

i 

0:03:25 

6,102 

8.7164 

45.951 

"356 

i 

0:03:27 

6 ,1.62 

8.7262 

54.049 

392 

i 

0:03: 29 

6,221 

8.7357 

62.147 

'354 

i 

' 0:03:52 

6,906 

8.8401 

70.242 

406 

l 

0:04: 22 

7,799 

8.9617 

78.331 

390 

i 

0:04:23 

i 

7,899 

8.9656 

86.402 

“345 

i 

0:04:49 

8,603 

9.0599 

94.387 
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i 

APPENDIX A-7 


CONVERSION OF TIMES -TO -"FAILURE DATA TO CYCLES-TO-FAILURE - 
FOR 97,500 PSI ALTERNATING STRESS AND STRESS RATIO 0.70 


Spec. 

No. 

Machine 

No. 

Time to 
Failure 

Cycles to 
Failure 

L °g e 

Cycles 

Median 

Ranks 

254 

3 

0:07:17 

12,964 

9.4699 

3.778 

423 

1 

0:07:18 

13,038 

9.4756 

9.151 

408 

1 

0:07:29 

13,365 

9.5004 

14.581 

379 

1 

0:07:45 

13,842 

9.5355 

20.024 

367- 

1 

0:08:20 

14,883 

9.6080 

25.471 

370 

1 

0:09:21 

16,699 

9.7231 

.30.921 

314 

3 

0:10:44 

19,105 

9.8577 

36.371 

237 

2 

0:11:19 

20,189 

9.9129 

41.823 

307 

3 

0:11:21 

20,203 

9.9136 

47.274 

, 292 

2 

0:11:28 

20,457 

9.9261 

52.726 

261 

2 

0:11:48 

21,004 

9.9526 

58.177 

256 

3 

0:12:59 

23,110 

10.0480 

63.629 

274 

3 

0:13:41 

24,356 

10.1005 

69.079 

262 

3 

0:13:51 

24,653 

10.1127 

74.529 

312 

3 

0:14:03 

25,009 

10.1270 

79.976 

258 

3 

6:14:08 

25,157 

10.1329 

85.419 

227 

3 

0:15:52 

28,243 

10.2486 

90.849 

295 

3 

0:18:03 

32,129 

10.3775 

96.222 


APPENDIX A-8 

CONVERSION OF TIMES-TO-FAILURE DATA TO- CYCLES -TO -FAILURE 


FOR 

76*500 PSI 

ALTERNATING 

STRESS AND 

STRESS RATIO 

0.70 

Spec. 

Machine 

Time to 

Cycles to 

Lo g e 

Median 

No. 

No. 

Failure 

.Failure 

Cycles 

Ranks 

382 

1 

0:22:09 

39,560 

10.5856 

3.778 

441 

1 

0:24:36 

43,936 

10.6905 

9.151 

416 

1 

0:26:01 

46,466 

10.7465 

14.581 

405 

1 

0:26:38 

47,567 

10.7699 

20.024 

471 

1 

0:29:30 

52,687 

10.8721 

25.471 

363 

1 

0:32:41 

58,372 

10.9746 

30.921 

141 

1 

0:33:35 

59,980 

11.0018 

36.371 

427 

1 

0:33:35 

59,980 

11.0018 

41.823 

285 

3 

0:34:43 

61,796 

11. ,0316 

47.274 

404 

. 1 

0:34:51 

62,242 

11.0388 

52.726 

200 

1 

0:36:52 

65,844 

11.0981 

58.177 

154 

1 

0:37:51 

67,491 

11.1214 

63.629 

277 

3 

0:37:55 

67,600 

11.1298 

69.079 

185 

1 

0:38:43 

69,184 

11.1440 

74.529 

218 

1 

0:39:18 

60,190 

11.1590 

79.976 

176 

1 

0:40:19 

72,006 

11.1845 

85.419 

326 

3 

0:42:35 

75,798 

11.2358 

90.849 

172 

1 

0:43:35 

77,840 

11.2624 

96.222 



APPENDIX A-9 


CONVERSION OF TIMES -TO -FAILURE DATA TO CYCLES -TO -FAILURE 


FOR 

70,000 PSI 

ALTERNATING 

STRESS AND 

STRESS RATIO 

0.70 

Spec • 

Machine 

Time to 

Cycles to 

Loge 

Median 

No. 

No. 

Failure 

Failure 

Cycles 

Ranks 

272 

3 

0:55:50 

' 99,383 

11.5067 

3.778 

385 

1 

0:56:29 

100,879 

11.5217 

9.151 

321 

3 

0:56:51 

101,193 

11.5248 

14 . 581 

360 

1 

0:58:09 

103,856 

11.5508 

20.024 

327 

.3 

0:58:38 

104,367 

11.5557 

25.471 

442 

1 

1:02:04 

110,851 

11.6160 

30.921 

268 

' 3 

1:06:16 

117,955 

11.6781 

36.371 

229 

3 

1:07:44 

120,565 

11.7000 

41.823 

304 

3 

1:12:36 

129,228 

11.7693 

47.274 

230 

3 

1:14:08 

131,957 

11.7902 

52.726 

231 

3 

1:16:23 

135,962 

11.8201 

58.177 

228 

.3 

1:17:09 

137,327 

11.8301 

63.629 

337 

1 

1:18:35 

140,350 

11.8519 

69.079 

413 

1 

1:19:44 

142,404 

11.8664 

74.529 

247 

3 

1:20:27 

143,201 

11.8720 

79.976 

287 

2 

1:24:03 

149,945 

11.9180 

85.419 

255 

3 

1:29:29 

159,280 

11.9784 

90.849 

232 

3 

• 1:34:14 

167,735 

12,0301 

96.222 



